Quantum search with trapped ions 
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We propose an ion trap implementation of Grover's quantum search algorithm for an unstructured 
database of arbitrary length N. The experimental implementation is appealingly simple because 
the linear ion trap allows for a straightforward construction, in a single interaction step and without 
a multitude of Hadamard transforms, of the reflection operator, which is the engine of the Grover 
algorithm. Consequently, a dramatic reduction in the number of the required physical steps takes 
place, to just 0(s/N), the same as the number of the mathematical steps. The proposed setup allows 
for demonstration of both the original (probabilistic) Grover search and its deterministic variation, 
and is remarkably robust to imperfections in the register initialization. 

PACS numbers: 03.67.Lx, 03.67.Ac, 03.67.Bg, 42.50.Dv 



Introduction. Quantum computers hold the promise 
of a dramatic speed-up for certain types of problems, 
which are processed intrinsically slowly by classical com- 
puters The fast search invented by Grover is ar- 
guably the quantum algorithm with the most outreach- 
ing implications. The Grover algorithm searches an ar- 
bitrary element in an unsorted database with N entries 
quadratically faster than its classical counterpart, with 
only 0(y/N) calls to an oracle. It is realized with the 
repeated application of two operations on an initial uni- 
form superposition of all states - an oracle, which flips 
the phase of the marked element \ip m ), and a reflection 
of the state vector about the mean. As N increases, the 
fidelity approaches unity, with error 0(1/N). 

Grover's algorithm has been demonstrated experimen- 
tally with two (N = 4) 3] and three (N = 8) QJ qubits 
in nuclear magnetic resonance, and with N = 32 items in 
classical optics [5] . This latter experiment demonstrated 
that entanglement is not essential for Grover's search but 
only the superpositional nature of the data register [f| . 

Trapped ions present one of the most plausible physical 
implementations of a quantum computer for their quan- 
tum dynamics is readily controlled and characterized and 
the system is scalable in principle. The theory 
number of landmark trapped-ion experiments [8j, 
been a major driving force behind quantum computing in 
recent years. Grover's search, however, has been imple- 
mented with a pair of ion qubits (N = 4) only jT^], largely 
because of the complexity of an ion-trap search experi- 
ment stemming from the necessity to construct the re- 
flection operator, the standard implementation of which 
requires numerous quantum (mainly Hadamard) gates. 

The reflection operation used in Grover's algorithm is 
known in matrix theory as Householder reflection (HR) 
and is widely used in classical data analysis. We have re- 
cently advocated the use of HR as a very efficient tool for 
a range of quantum-state engineering problems in single- 
particle [HI, HH and many-particle systems [HI ■ These 
include the synthesis of arbitrary unitary transformations 




superposition-to-superposition transitions [l3 ]. and 
engineering of arbitrary mixed states [l2j | of a single par- 
ticle, and quantum-state engineering of collective entan- 
gled states of trapped ions [13| . This work has been trig- 
gered by the fact that the HR operator emerges naturally 
as the propagator in a coherently-driven degenerate two- 
level system: as a single HR for an iV-pod sy_stem of one 
degenerate and one non-degenerate level [14], and as a 
product of HRs in a system of two degenerate levels [l5| . 

In this Letter, we make use of the HR implementa- 
tion with collective states of trapped ions [131 ] to propose 
a very simple experimental demonstration of Grover's 
search. The linear ion-trap architecture and the ensuing 
linkage pattern prove ideally suitable for Grover's search, 
because the required HR operator arises naturally, as a 
solution of the Schrodinger equation for a Hamiltonian 
inherent for laser-driven trapped ions. Hence the global 
HR operator, used to perform the reflection about the 
mean, can be produced in a single interaction step, rather 
than by O(N) gates hitherto, which greatly reduces the 
required overall nonquery resources [1 61 ] . 

As stressed by Grover 16], the total number (query 
and nonquery) of physical steps in implementations of the 
reflection operator with a multitude of Hadamard gates 
can exceed considerably the number of logical steps. For 
example, in one of the most sophisticated experiments 
so far [J], about 100 pulses have been used to simulate 
Grover's search in an TV — 8 register. The implementa- 
tion proposed here is the first one, where the total num- 
ber of physical steps is the same as the number of algo- 
rithmic steps, which is N G = [n/ (2 sin" 1 (2y / N - 1/N))] 
([n] denotes the integer part of n), or Nq — [(tt/4)Vn] 
for large N. This speed-up reduces significantly the ex- 
periment complexity, the overall processing time, and the 
fidelity requirement for each interaction step. 

Hamiltonian. We consider a string of N two-state 
ions confined in a linear Paul trap. The radial trap fre- 
quencies are much larger than the axial frequency w z , so 
that the ions form a linear string along the z axis [l7j |. 
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tion the Hamiltonian |l| can be written as [H 



FIG. 1: Linkage pattern of the collective states of N trapped 
ions driven by a red-sideband laser pulse. The lower state 
|V>o)|l) involves one vibrational phonon and no ionic excita- 
tions, whereas the upper set of states {|^>fc)|0)}*^=i involves 
zero phonons and a single ionic excitation of the kih ion. 



The ions are laser-cooled into their motional zero-phonon 
ground state |0) [l8| and their internal dynamics is con- 
trolled by a laser (or two lasers for a Raman transition), 
with a carrier frequency tuned near the red-sideband res- 
onance: Lu n = luq— uj z +5 n , where loq is the Bohr transition 
frequency of the ions and S n is the detuning of the nth ion 
(\S n \ <C ojz)- In the Lamb-Dicke limit and the rotating- 
wave approximation, the interaction Hamiltonian is 17 1 
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(1) 



Here and o are the creation and annihilation opera- 
tors of the center-of-mass phonons, while er+ = |l n )(0 n | 
and cr~ = |0„)(1„| are the Pauli spin-flip operators 
for the internal states |0„) and |1„) of the nth ion. 
rj n = x/frk^ cos 2 6 n /2Moj z is the Lamb-Dicke parame- 
ter, where k n is the laser wavevector, 9 n is the angle be- 
tween the trap axis and the laser propagation direction, 
M is the ion mass, and f2„ (t) is the time-dependent real- 
valued Rabi frequency of the laser-ion coupling. All Rabi 
frequencies must have the same time-dependent envelope 
/ (i), but they may have different amplitudes f2„. 

Since the Hamiltonian ([T]) is of Jaynes-Cummings type, 
it conserves the sum of ionic and vibronic excitations; 
hence the Hilbert space decomposes into subspaces with 
a definite number of excitations. The single- excitation 
subspace is spanned by the N + 1 basis states, 



|^)|0) = |0l 
IV-o>|1) = |0i- 
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where is the vibrational state with /i phonons (u = 
0, 1), and |0i . . . 0fc_i, 1&, Ok+i ■ ■ - On) is a collective ionic 
state (fc = 1, 2, . . . TV), in which the fcth ion is in state |1) 
and all other ions are in state |0). Our quantum memory 

register will be the set of states pa)) : {\i/>k)}k=ii through- 
out the algorithm operation the population remains in 
vibronic state |0); it will be omitted hereafter. 

Householder reflection. After a phase transforma- 



Hj(t) = ~ $>»(f)|^>ty»| + ~^><^| + h.C. 



(3) 
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where all detunings are assumed equal, S n = S, and 
9n(t) = r] n Q n e^' l ' n f(t)/y/N = g n f(t) is the coupling 
between the internal and motional degrees of freedom 
for each ion (n = 1, 2, ... , A). The linkage pattern for 
this Hamiltonian is an A-pod, wherein state |^>o) is cou- 
pled to each state in the manifold {\ipk)}^ = i, a s shown in 
Fig. [1] It has been shown very recently that the prop- 
agator within the single-excitation manifold {\tpk)}k=i 
for exact resonance (S = 0) and for root-mean-square 
(rms) pulse area A = g f(t)dt — 2(21 + l)n (with 



I = 0,1,2,.. 
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M(x) = l-2| X )(x|, 



is given by a 



(4) 



where the components of the A-dimensional normalized 
vector \x) are the normalized couplings, 



1 T 
\x) = ~ [9i,92,---,9n] 
9 



(5) 



The generalized HR within the manifold {\ipk)}k=f 

M(x;^) = l + (c^-l)| x )(x|, (6) 

is realized at suitably chosen detunings. For example, for 

a hyperbolic-secant pulse, f(t) = sech(i/T), with rms 

area A = 2-nl (1 = 1,2,...), the HR phase ip is produced 

l-l 

by detunings, obeying the equation 2 arg Y\ [ST + i(2j + 



1)] = tp [14{. The generalized HR can be created also for 
other pulse shapes, e.g. Gaussian, but the required pulse 
area and detuning have to be evaluated numerically. 

Grover search. The implementation of the Grover 
algorithm is depicted in Fig. [21 

1. Initialization We first prepare the register (I2a|) in 
an equally weighted superposition of all A states |2J, i.e., 
\i>w) = £n=i IVO/VA =[1,1,..., 1] T /VN; this is the 
well-known entangled W-state, which has been demon- 
strated experimentally with eight ions For this, one 
can choose between the latter technique adiabatic- 
passage techniques IjJ 2(| , or a HR-based technique 13 1 . 
For Grover's search it is most convenient to use the fol- 
lowing alternative [l3 |: prepare the ion chain initially 
in state |^o)|l)i an d then apply, uniformly at all ions, a 
red-sideband pulse with a temporal area of tt/vA expe- 
rienced by each ion (hence rms area of 7r) . 

2. Grover iteration Each Grover logical step consists of 
two operators, an oracle call and a global HR 2]. Firstly, 
the oracle marks the searched state \ip m ) by inverting its 
phase, which in our implementation is achieved by indi- 
vidual addressing of only the mth ion by a local 2it pulse; 
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FIG. 2: (Color online) Implementation of Grover's search with 
N trapped ions. The ion chain is initialized in an entangled 
IF-state by preparing it initially in state |i/)o)|l), and then 
apply a red-sideband pulse with a temporal area of n/\fN 
uniformly at all ions. Then we perform repeatedly, Ng times, 
a local reflection of the marked ion by applying a 2n pulse on 
it, followed by a global reflection, inflicted by a 2ir/y~N pulse 
applied to the entire ion string. Finally, the marked ion is 
identified for it is the only one in state 1 1) . 



in fact, this is a HR operation M(% m ) with an interaction 
vector \xm) identical to the basis vector \ip m ) of Eq. (|2ap. 
The second HR M.(xw) with the HR vector \xw) = \^w) 
inverts the amplitude of all states \ipk) about the mean. 
In our implementation this HR is produced by applying 
a single laser pulse, uniformly at all TV ions in the regis- 
ter, with a temporal area of 2ir/vN (which amounts to 
rms-area of 2tt), as shown in Fig. [2] After the execution 
of the Grover operator Nq times, the system is driven 
into the marked state, [M(x^)M(x m )] G \tpw) ~ \i>m)- 

3. Detection The marked ion is identified for it is the 
only one in state |1). 

In Fig. [3] the probability to find the searched state for 
N = 15 ions is plotted as a function of time. In three 
interation steps, which amount to 3 local (oracle) pulses 
and 3 global (reflection) pulses, this probability increases 
to about 0.92, after which it decreases (as a part of os- 
cillations between zero and unity in a long run). 

Deterministic Grover search. The original Grover 
algorithm is probabilistic (except for N = 4): it finds 
the marked state with a probability (fidelity) close to, 
but less than unity. A slight modification of the origi- 
nal Grover algorithm, with the supplement of the (real) 
reflection operators by complex phase factors, makes the 
search fully deterministic, with a unit fidelity for any TV, 
with Nq or N G + 1 roundtrips [2l| . 

The deterministic search can be implemented with 
the same strategy as the original Grover search de- 
scribed above, by replacing the HRs M(x m ) and M.(xw) 
with the generalized HRs M(x m ;<p) and M(xn/;^): 

[M(xw;v)M( Xm ]V)] Na \ifar) = \i>m)- The HR phase 
reads ip = 2 sin -1 [ViVsin(7r/(4/V G + 6))] [21]; it can be 
produced by a suitable detuning 8, as explained above. 
In Fig. [3] the population of the marked state is plot- 
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FIG. 3: (Color online) Numerically calculated population of 
the marked state vs time for N = 15 ions for probabilis- 
tic (dashed line) and deterministic (solid line) Grover search. 
The upper frame depicts the sequence of local (dashed, Ok, 
area 2n) and global (solid, I k , area 2tt/^/N) pulses, both with 
sech shapes. The individual couplings g n are given by the 
components of the QHR vectors \xm) and |xw), each multi- 
plied by g = 2/T. The detuning in the deterministic search is 
ST ~ 0.589, which produces the required phase ip ~ 0.6617T. 
Maximum probabilities occur after the completion of I3, when 
the dashed curve reaches 0.92 and the solid line unity. 



ted for deterministic Grover search with N — 15. After 
just three iterations, the occupation probability for the 
marked state approaches unity, whereas it is only 0.92 
for the standard Grover search. 

Practical considerations. Now we discuss briefly 
various issues that may arise in a real experiment. An ob- 
vious deviation from the idealized theory can be the im- 
perfect initialization of the search register (the W-state) , 
due to an inhomogeneous spatial profile of the driving 
laser(s), for example, if the laser beam for the global re- 
flection is tilted from being collinear with the trap axis. 
Then the initial state will not be an equal superposition, 
but the more general state |\I/ a ) = Yln=i ^li^n), where 
the probability amplitudes a n may deviate from l/s/N, 
and they may also be complex, due to possible phase 
differences between the ions. It has been shown that 
the register |^ ) can still be used for Grover's search if 
the amplitudes a n do not deviate too much from l/\fN 
[22I . l23| . Biham et al. (23[ have shown that the appli- 
cation of the same Grover operator, as for an equally 
weighted initial superposition, to the initial state |\l/ a ) 
still produces the marked item with a reasonable proba- 
bility. However, it is readily shown that a HR with the 
same vector as the amplitude distribution in \^ a ), 



\x) = [ai,a 2 , ...,a N ] 



(7) 



allows for a higher probability. Using the approach of 
Hsieh and Li [2 41 ] one can show that the optimum number 
of steps Nq is replaced by N a = [7r/(4|a m |)] (for small 
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FIG. 4: (Color online) Infidelity of finding ions 1, 5 and 10 
after 3 steps for N = 20 for a Gaussian spatial laser beam 
profile of the global pulse (rms area 2tt, detuning 8 = 0), 
where the outside ions experience only a fraction 1 — s of the 
laser intensity in the middle of the ion chain. 



a m ), where a m is the initial amplitude of the marked 
state \ipm). Obviously, for N a to apply to any state in 
the register \^ a }, the largest and smallest amplitudes in 
Eq. ([7]) should produce approximately the same N a . 

Figure 0] shows the infidelity to detect the marked 
states | -0i ) , \ips) and |V>io) for N — 20 ions after Nq = 3 
steps for a Gaussian laser beam profile of the global pulse, 
where the outside ions experience only a fraction 1 — e of 
the laser intensity in the middle of the ion string. The 
proposed technique proves to be very robust with respect 
to the imperfectness of the laser profile. 

In the proposed implementation the interaction vec- 
tor |(7|) can be produced very simply: by using the same 
laser beam as in the initialization step (which has pro- 
duced the possible inhomogenuities in the register), but 
with twice as high Rabi frequency, cf. Fig. [51 This ro- 
bustness to imperfections in the initial register, caused by 
unequal individual laser-ion couplings, also implies that 
the proposed implementation of Grover's search can uti- 
lize higher vibrational modes, e.g., the breathing mode, 
for which the laser-ion coup lings depend on the position 
of the ion in the string [171 ] ; using higher phonon modes 
greatly reduces deleterous heating effects. 

Discussion and conclusions. In this Letter, we 
have proposed a very simple and concise implementa- 
tion of Grover's search algorithm in a linear chain of N 
trapped ions. Unlike earlier proposals, which, in addition 
to the 0(y/N) queries, require many more physical non- 
query steps, the proposed implementation requires only 
0(y/N) such steps, one for each of the 0(yN) queries; 
hence each Grover iteration is performed by one local 
HR (query) and one global HR (nonquery). The speed- 
up in regard to physical steps derive from the natural 
emergence of the HR operator as the propagator in a 



red-sideband laser-driven linear string of trapped ions. 

We point out that the proposed implementation, while 
using an entangled W-state as a quantum register, does 
not utilize the full power of a quantum register, which can 
contain 2 N states. Work in this direction, which requires 
dealing with products of HRs, is in progress [251 ] . 
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the Bulgarian National Science Fund under grants VU- 
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